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Введение. Сопоставлены два фактора токарной обработки 
металлов резанием: рассеиваемая температура и мощность 
необратимых преобразований в материале изделия и ин- 
струмента. Статья посвящена вопросам математического 
моделирования их связи. 

Материалы и методы. Математический аппарат основан 
на модификации уравнения Вольтерры, которая предпола- 
гает использование двукратного интеграла. Он показывает, 
как выделенная ранее при резании тепловая энергия влияет 
на текущее состояние температуры в зоне контакта ин- 
струмента с обрабатываемой деталью. В работе помимо 
предложенной новой базовой математической модели ис- 
пользованы результаты обработки экспериментальных 
данных о мощности необратимых преобразований и об 
измеряемой температуре в зоне контакта инструмента с 
деталью при токарной обработке металла. Эксперименты 
проводились на станке 1К625 и стенде ЗТО.201-1. В пакете 
МаЧаб была создана специализированная программа для 
обработки массивов информации, описывающих процессы, 
протекающие при резании (силы реакции, вибрации ин- 
струмента и мощность необратимых преобразований). В 
этой же программе выполнен расчет температуры в зоне 
контакта инструмента и детали. 

Результаты исследования. Проведена процедура парамет- 
рической идентификации предлагаемой базовой математи- 
ческой модели. Полученная модель показала высокую сте- 
пень близости экспериментальных данных о состоянии 
температуры в зоне резания и моделируемого уровня теп- 
ловой энергии, однако на начальном участке измеряемой 
температурной зависимости результаты двух этих подхо- 
дов не вполне согласуются. Это можно объяснить погреш- 
ностью экспериментального измерения температуры на 
основе оценки вырабатываемой термоэлектродвижущей 
силы (термоЭДС), которая формируется в результате обра- 
зования естественной термопары в зоне контакта инстру- 
мента и детали. 
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Обсуждение и заключения. Предлагаемая математическая 
модель позволяет адекватно описывать процесс преобразо- 
вания механической составляющей энергии резания в теп- 
ловую составляющую через показатель суммарной выделя- 
емой мощности механического взаимодействия в зоне ре- 
зания за все время обработки. 


Ключевые слова: точение, резание, мощность необрати- 
мых преобразований, обработка металлов, теплопровод- 
ность металлов, уравнение Вольтерры, термоэлектродви- 
жущая сила (термоЭДОС). 


Образец для цитирования: Бордачев, Е. Математическое 
моделирование температуры в зоне контакта инструмента и 
изделия при токарной обработке металлов / Е. Бордачев, 
В. П. Лапшин // Вестник Дон. гос. техн. ун-та. — 2019. — 


Б5си5яоп апа Сопсиз1оп5. ТВе ргорозе4 тафетайса| по4е! 
епа ез ю адедиае1у 4езстфе е сопуег1оп оЁ Фе теспашса| 
сотропеп{ ог фе сийш» епеггу шю фе Феппа! сотропеш 
@гоцэВ фе ш@саюг оЁ Фе 1о{а| ошрий ро\ег оЁ Фе тесвап- 


са] ифегасноп ш Фе сито топе Юг аП Ше ргосеззте те. 


Кеуумог@$: шготе, сибе, роууег оЁ итеуегзе ЧтапзЮт- 
таНопз, тейа| ууогкте, Фегта[ сопдаспуйЙу оЁ теа15, УоКег- 
га ищерта| едиаНоп, КепповесйотоНуе гсе (фегта! ЕМЕ). 


Еог сйайоп: Е. Вогдавеу, У.Р. ГарзЬш. Мафетанса| {ет- 
регабаге зииайоп ш 1ю0]-ю-уогК сошас{ хопе файпе шей] 
баготе. Уезии К оРГОЗТО, 2019, уо1. 19, по. 2, рр. 130-137. 
Ьрз://401.0г2/10.23947/1992-5980-2019-19-2-130-137 


Т. 19, №2. — С. 130-137.  Б@рвз://401.ог/10.23947/1992- 
5980-2019-19-2-130-137 


Введение. Тепловые процессы, протекающие при обработке металлов на металлорежущих станках (как 
в металлах, так и в других средах), являются сложноописываемыми. Их анализ основан на уравнении тепло- 
проводности Фурье [1]. В случаях, связанных с обработкой металлов резанием, нужно иметь представление о 
температуре в зоне контакта инструмента и обрабатываемой детали [2—5]. При этом нет необходимости кон- 
тролировать распространение тепла за пределами данной зоны, так как на условиях обработки это не сказыва- 
ется. Аналогичные рассуждения можно применить к другому широко распространенному методу описания 
тепловых процессов, который учитывает формирование истоков и стоков тепла в детали [6]. Такие подходы не 
позволяют использовать динамику процесса резания в описании возникающих тепловых явлений в зонах кон- 
такта и создать систему уравнений, представляющих преобразование механической энергии в тепловую. 

При операциях с металлами на металлорежущих станках необратимые преобразования в детали и ин- 
струменте вызывают производство тепла, которое, согласно второму закону термодинамики, рассеивается в 
пространстве. Тепловые процессы в металлах более инерционны, чем динамика процесса обработки, поэтому 
рассеивание продолжается достаточно длительное время. Само производство тепла в рассматриваемом случае 
связано с мощностью необратимых преобразований [11], которая зависит от величины сил реакции на формо- 
образующие движения инструмента [2—5] в зоне контакта с деталью. При этом производство тепла является 
основным механизмом рассеивания выделенной энергии преобразования материалов [12]. 

Динамика процесса резания и вопросы ее моделирования напрямую связаны с представлением сил и 
реакций в координатах состояния процесса. В этом случае, как правило, учитывается только механическая со- 
ставляющая, что видно из зарубежных [13—18] и отечественных работ [19—21]. Таким образом, целесообразно 
ввести в моделирование динамики процессов обработки металлов координаты, описывающие изменение тем- 
пературы в зоне контакта инструмента и детали через учет необратимых преобразований. Такой подход позво- 
лит качественно повысить адекватность используемого математического аппарата. 

Материалы и методы 

Обоснование математической модели. Силы, препятствующие внедрению инструмента в материал 
обрабатываемой детали, имеют сложное пространственное расположение в системе координат, связанной с 
осями деформации инструмента. В литературе, посвященной анализу сил реакций, возникающих в зоне реза- 
ния, принято разложение такой реакции на составляющие [4—6, 10, 19—21]. В рассматриваемом нами случае 
будут справедливы следующие соотношения между составляющими сил реакции [5]: 


2 2 2 
В=Р2+Р2+Р2, (1) 


где в — общая силовая реакция со стороны процесса резания на формообразующие движения инструмента; Р. 
—Щ проекция силовой реакции на ось главного вращательного движения резания в направлении подачи; Р, === 


проекция силовой реакции, направленная по радиусу главного вращательного движения резания в вершине ре- 
зания, В, —Й проекция силовой реакции на ось, совпадающую со скоростью главного движения в вершине ин- 


струмента. 
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Соотношения между введенными силами Р,,Е,,Ё, зависят от многих факторов: геометрия резца, сте- 


пень износа резца и т. д. [4]. Так, в [5] при обработке острым резцом с параметрами у = 15°, ф= 45° и/,= 0° со- 
отношение в среднем равно 
РРР. = (0,3 — 0,4), (0,4 — 0,5). (2) 
Аналогичные рассуждения можно провести по анализу скоростей относительного движения инстру- 
мента вдоль обрабатываемой детали и получить общую скорость такого движения: 


И. +И2+И.2. (3) 


С учетом этих рассуждений определим мощность необратимых преобразований: 


[ы 8 2 2 [/2 2112 
№М=уЬ/ +В +Е/ 7; +, +Т,. (4) 


Допустим, что производство тепла при резании в металлорежущих станках линейно зависит от значе- 
НИЯ МОЩНОСТИ необратимых преобразований: 
Т, = М. (5) 


Здесь ТГ, — величина, описывающая приращения температуры в зоне резания (°С); Кг — коэффициент связи 


между мощностью необратимых преобразований в зоне резания и величиной приращения тепловой энергии в 
зоне контакта инструмента и детали размерностью с°С/Нм. Однако общее значение температуры в заданный 
момент времени и в заданной зоне контакта инструмента и детали будет определяться не только текущим зна- 
чением прироста тепла. В связи с рассеиванием тепла в материале детали и инструмента следует учитывать 
влияние: 

— мощности необратимых преобразований по всему предыдущему пути, 

— времени обработки изделия. 

Представим себе обработку в виде некоторого дискретного графика процесса, где (и) — текущее зна- 
чение мощности необратимых преобразований; Г[(п) текущее значение пути, пройденного инструментом при 
обработке; #(п) — текущее значение времени обработки. Согласно (5), в каждой рассматриваемой точке дис- 
кретного графика процесса обработки будет выделено приращение тепла, пропорциональное мощности необра- 
тимых преобразований в этой точке пространства и времени. Текущее значение температуры в точке и будет 
определяться как приростом тепла в этой точке, так и влиянием приростов температуры, происходивших ранее 
на пути, пройденном инструментом (Г) во время обработки (1). С увеличением времени и пройденного пути 
ослабевает влияние этих факторов на процесс выработки тепла в рассматриваемой зоне контакта, поэтому для 
математического описания удобен подход, основанный на использовании оператора Вольтерры второго рода 
[10]. Его невозможно применить напрямую из-за сложности описания распространения тепла в металлах, по- 
этому в качестве базовой модели примем мультипликативный критерий оценки влияния предыдущего прироста 
тепла на текущее значение в виде двойного интеграла: 

7, = ©, + Аг [[,@-Ви(п-ОмМтачам. (6) 

р 
Здесь Т, — значение тепла в зоне контакта инструмента и обрабатываемой детали; ®; — температура окру- 
жающей среды; иг (у-Г) — ядро оператора, характеризующее влияние выделенной ранее мощности необра- 
тимых преобразований по пути обработки на текущее значение температуры; и",(п-#) — ядро оператора, ха- 
рактеризующее влияние выделенной ранее мощности необратимых преобразований по времени обработки на 


текущее значение температуры; р — область интегрирования, характеризующая пространственно-временное 
рассеивание тепла; ‘у имеет размерность расстояния, м; | — имеет размерность времени процесса обработки, 


с. Сами ядра интегрального оператора представляют собой безразмерные коэффициенты усиления. 
Интегральный оператор, представленный в виде двойного интеграла (6), можно свести к кратному ин- 
тегралу следующего вида: 
(и) 1 


7. -0; + | ма-В ау [м (т-0мая. (7) 
0 0 


Пройденный инструментом путь является функцией времени, т. к. каждому значению времени можно 
поставить в соответствие значение пройденного пути. Следовательно, подход (7) справедлив. 
В качестве ядер предложенного в (7) интегрального оператора удобно использовать экспоненциальные 
функции следующего вида [10]: 
р-р) =е* 2 
и, (п-1) = г’ 


где 01, 06 — масштабирующие параметры интегрального оператора, подлежащие идентификации. 


(8) 
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Таким образом, предложенный в (7) интегральный оператор примет следующий вид: 
770) 1 
Т, 9,4 | с“ Эду ЭмМтат. (9) 
0 0 
Этот оператор имеет решение для стационарного случая, когда мощность необратимых преобразований 
является постоянной величиной № = М(®. 


Ат № 
0195 


т. =©, + (-е “4-е. (10) 


Как видно из (10), при стационарном процессе обработки температура в зоне резания экспоненциально 
нарастает с начального значения ©, до некоторого конечного значения, равного: 


Ат. № 
0192 | 


Т. =©, + (11) 


Проиллюстрируем эти рассуждения, задав ®, =25 °С; Аг= 0,0026 с°С/Нм; ви =0,03; @> =0,01; 
№ =20 Нм/с (рис. 1). 


Т2(С°) 


Им) 





Рис. 1. Реализация процесса стационарного нарастания температуры в зоне резания 


Как видно из рис. 1, стационарный тепловой процесс в зоне резания нарастает экспоненциально и, до- 
стигнув некоторого равновесного состояния, больше не изменяется. Реальные процессы, протекающие в ме- 
таллорежущих станках, не являются стационарными [11-12]. Исходя из этих соображений, процесс, представ- 
ленный на рис. |, можно рассматривать как движение, к которому в устойчивом случае стремится температура 
зоны контакта инструмента и детали. При этом уравнение (10) является порождающим уравнением для инте- 
грального оператора (9). 

Результаты исследования 

Параметрическая идентификация модели по экспериментальным данным и итоги моделирова- 
ния. Для параметрической идентификации модели проведен натурный эксперимент на токарном станке 1К625 
с установленным стендом 5ТО.201-1. Стенд предназначен для изучения динамических и тепловых процессов, 
протекающих при резании металлов на токарных станках в различных режимах. 

Основной целью эксперимента было исследование влияния динамики резания на тепловые процессы в 
зоне контакта инструмента и детали. Под изменением динамики резания понимается взаимное влияние износа 
инструмента и процесса резания, связанное с возникающей в зоне резания силовой реакцией. Измерение этой 
реакции и температуры через термоэлектродвижущую силу (термоЭДС), формируемую в зоне контакта ин- 
струмента и детали, позволяет оценить влияние динамики резания на тепловые процессы в указанной зоне. 

Во время эксперимента выполнялось непрерывное точение детали типа вал. При этом сохранялись все 
параметры обработки и осуществлялся контроль изменения силовой реакции со стороны процесса резания. 
Условия и основные результаты эксперимента подробно описаны в [22]. 
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Данные, полученные в результате натурного эксперимента, носят дискретный характер, следовательно, 
к ним нельзя напрямую применить модель, представленную интегральным оператором (9). По этой причине 
приходится видоизменить базовую математическую модель, описывающую динамку тепловых процессов в 
контакте инструмента и детали. Для начала оставим в рассуждениях только член интегрального оператора, за- 
висящий от мощности необратимых преобразований: 


Ци) 1 
ТА =Ю | ау е 0 М(мат. (12) 
0 0 


Так как эксперимент уже проведен, то время его окончания | и соответствующее этому времени зна- 
чение пути Ца.) являются постоянными величинами. Следовательно, будет справедливо вынести экспоненты 


с этими значениями за пределы интегралов: 
ЩЕ) | 
ТМ = те ео, | ета | е- 1 М(т)ат . (13) 
0 0 
Примем, что существует некоторое дискретное множество расчетных значений мощности необратимых 
преобразований, полученных в результате анализа дискретной экспериментальной информации о динамике 


процесса резания, то есть №(#) = М, М, №, №4......... Исходя из этого, интегральный оператор (13) примет вид 
следующей суммы: 
Г. 1 Г. я 
ТМ = те ей [м | су | е°-Пам + № | су | °-Пат+ > 
2. й ° . й 1, . - (14) 
> +М; | е* у | с°-Пат+ №4 | ету [ с°Пат+ вы 1] 
Г, 1, Г, ь 


Сделав допущение ЛМ; =0, получим окончательное решение для случая, описываемого выражением 
(14): 
ти - АВА [5 (с) — 9) (ебыь — еб + «> 
о: (15) 
<> +М№;3 (Е) — еб.) еб еб)... ] 
Выражение (15) очень удобно для обработки больших массивов дискретной информации о процессе 
резания. Фактически речь идет о взвешенной относительно конечного значения сумме мощностей необратимых 





преобразований. 

Рассмотрим только один эксперимент из серии, при котором износ по задней поверхности инструмента 
отсутствовал. В следующих случаях менялась изношенность инструмента, а значит, и динамика процесса (как с 
механической, так и с термодинамической точки зрения). Отметим, что полученные в этом эксперименте ре- 
зультаты обобщаются для всех остальных случаев, однако объем статьи не позволяет их рассмотреть. 

Реакция со стороны процесса на формообразующие движения инструмента вдоль оси х представлена 
на рис. 2. 
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Рис. 2. Реакция со стороны процесса на формообразующие движения инструмента вдоль оси х на первом шаге эксперимента 
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Как видно из рис. 2, реакция на формообразующие движения инструмента при токарной обработке но- 
сит существенно нестационарный характер, но при этом наблюдается определенное притягивающее многообра- 
зие, вокруг которого и происходит возмущенное движение инструмента. Аналогично представленному на 
рис. 2 случаю в эксперименте получены реакции на формообразующие движения инструмента вдоль остальных 
осей. Кроме того, доступна информация об ускорениях инструмента относительно обрабатываемого изделия. 
Интегрирование этих данных с учетом известных элементов резания позволило получить значения скоростей 
обработки и пути, пройденного инструментом вдоль обрабатываемой детали. 

В пакете МаНаб была создана специализированная программа для обработки массивов информации, 
описывающих процессы, протекающие при резании (силы реакции, вибрации инструмента и мощность необра- 
тимых преобразований). В этой же программе выполнен расчет температуры в зоне контакта инструмента и 
детали по (15). Параметры модели представлены в табл. 1. 








Таблица 1 
Параметры идентифицированной модели 
1 02 Ат с°С/Нм ©, °С 
0,00099 0.00078 0.000159 24,8 




















На рис. 3 представлены результаты моделирования выражения (15) с учетом параметров табл. 1, а так- 
же экспериментально измеренная характеристика. 
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Рис. 3. Сравнение результатов, полученных в эксперименте и рассчитанных по модели 


Сопоставление рис. 3 и рис. 2 показывает, что измеренное значение температуры при врезании инстру- 
мента в деталь увеличивается с 25 °С (температура окружающей среды) до 330-340 °С. Согласно данным рис. 3, 
сначала наблюдается более быстрый рост измеряемого значения термоЭДС, а в конце процесса термоЭДС ста- 
билизируется. Это объясняется тем, что на стенде 5ТО.201-1 измеряется не сама температура в зоне контакта, а 
термоЭДС естественной термопары, образуемой в контакте инструмента и детали. В начале эксперимента из-за 
резкого скачка температурных градиентов термоЭДС показывает слишком большой и быстрый рост темпера- 
туры, а затем, по мере расхождения температурных градиентов стабилизируется. Этот эффект широко известен. 
Например, в [5] сказано: «Метод естественной термопары ненадежен, так как температура резания стабилизи- 
руется в течение 2-3 с». 

С точки зрения динамики, изменение температуры (рис. 3), полученное в результате моделирования 
выражения (15), близко к процессам, протекающим при резании (рис. 2). На графике виден продолжающийся 
рост температуры, обусловленный продолжающимся ростом сил реакции (рис. 2). Также график отражает все 
динамические особенности процесса, в том числе малоамплитудные скачки сил. Из этого можно сделать сле- 
дующий вывод: предлагаемая модель более точно описывает изменения температуры резания, чем данные 
естественной термопары. Также из рис. 3 ясно, что моделируемая характеристика, описывающая изменение 
температуры в зоне резания, в целом качественно совпадает с экспериментальной зависимостью. Эксперимен- 
тальная и моделируемая характеристики наиболее полно совпадают после стабилизации процесса резания, в 
пределах 15—24 секунд. 
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Обсуждение и заключения. Итоги исследования дают ответ на вопрос, каким образом можно добиться 
адекватного экспериментальным данным описания наследственного характера передачи и накопления тепла в 
процессе обработки металлов точением. С этой целью необходимо смоделировать температуру в зоне контакта 
инструмента и детали на основании предложенного модифицированного оператора Вольтерры. 

С научной точки зрения, наибольший интерес представляет предлагаемая базовая математическая мо- 
дель и результаты ее идентификации на основе эксперимента по реальной обработке металлов. Здесь наблюда- 
ется высокая степень совпадения моделируемых значений тепловой энергии в зоне контакта инструмента и об- 
рабатываемой детали в период установившегося процесса обработки. 
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